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ABSTRACT 


Theoretical  studies  have  been  carried  out  to  investigate  the  scattering  of  light 
by  spherical  aerosols  with  the  objective  to  obtain  basic  information  useful  for 
practical  applications.  Numerical  data  for  Mie  scattering  have  been  analyzed  to 
find  general  trends,  to  simplify  interpolation  problems,  and  to  establish  simple 
relationships  between  important  parameters.  Results  of  the  research  are  given 
in  the  form  of  abstracts  of  10  scientific  reports  and  6  related  papers  issued 
during  the  course  of  this  contract.  Revisions  and  additions  to  each  report  are 
given.  An  extensive  bibliography  of  tables  concerning  Mie  scattering  and  an  atlas 
of  scattering  diagrams  for  six  refractive  indices  from  n  =  1.  1  to  1.  5  and  a=  0.  5 
(0.  5)  10  are  given  as  appendixes. 
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Another  problem,  which  has  been  tackled,  concerns  the  absorption  and  scattering 
coefficients  for  small  aerosols.  Exact  derivations  of  the  equations  for  lde), 

K^a),  and  have  been  obtained.  The  asymptotic  expansion  of  the  exact  solu¬ 
tion  is  given,  where  the  leading  term  is  identical  with  the  Rayleigh  approximation. 


For  nonabsorbing  spheres,  the  total  scattering  coefficient  for  small  size  param¬ 
eters  is 


(5) 


For  all  practical  purposes,  the  error  remains  below  +_  1  percent  to  n  =  1.  4 
and  a  <1.  For  n-«  2,  the  error  reached  15  percent  at  a  =  1.  Figure  3  in 
Scientific  Report  3  shows  the  error  as  function  of  a  and  n  and  can  be  used  as 
a  reliable  guide  in  determining  the  upper  limits  of  a  and  n  for  which  equation 
(3)  is  useful. 


For  absorbing  spheres,  the  complex  index  of  refraction 

n  =  n  -  ix  (6) 

has  to  be  used,  and  the  total  extinction  coefficient  for  small  size  parameters 
is  given  by 


♦ 
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« 
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New  diagrams  of  the  angular  position  of  maxima  and  minima  have  been  con¬ 
structed  in  the  range  n  =  1.  05  to  1.  5,  where  the  data  for  n  =  1.  05  to  1.  3  are 
taken  from  Pangonis  and  Heller  (I960),  those  for  n  =  1.  486  from  Rowell 
(1963),  and  the  rest  from  the  Penndorf  and  Goldberg  tabulations  (1953). 

Rowell's  data  are  extremely  valuable  because  they  are  based  on  computations 
in  steps  of  Ad  =  1  degree.  However,  they  are  available  only  in  diagrams  and 
had  to  be  transformed  into  their  system.  For  ip  the  position  of  the  trenches 
could  also  be  determined  from  his  diagram.  Several  samples  are  selected  from 
the  present  collection  of  diagrams  of  angular  position  of  maxima  and  minima 
(figures  2  and  4)  and  modified  altitude  charts  (figures  3  and  5).  For  n  <  1.  33, 
both  the  maxima  (solid  line)  and  minima  (dashed  line)  are  shown;  for  larger 
refractive  indices  only  the  minima,  otherwise  the  diagrams  would  appear 
overloaded.  The  results  for  the  complete  collection  are  discussed,  but  it  is 
not  believed  necessary  to  show  them  all. 

For  i^,  a  <  7,  and  low  refractive  indices  (n<  1.  15),  only  diffraction  plays  a 
role;  and  the  system  of  maxima  and  minima  presents  itself  in  a  very  simple 
fashion  not  only  in  figure  3  but  also  in  figure  2.  For  n  =  1.2,  the  first  dis¬ 
turbance  in  the  diffractive  system  is  noted;  i.  e.  ,  the  effect  of  internal  reflec¬ 
tion  and  refraction  (reflective  system  for  short).  This  disturbance  begins  in 
the  backscattering  area  and  becomes  much  clearer  for  n  =  1.  25  and  1.  3  than 
for  n  =  1.  2,  and  finally  for  n  >1.  3,  only  the  first  diffractive  minimum  remains 
undistrubed. 

The  reflective  system  starts  for  such  refractive  indices  in  the  backward  area 
and  moves  toward  larger  scattering  angles  as  figure  2  indicates  and  more  or 
less  horizontal  in  the  altitude  charts  (figure  3).  The  case  for  n  =  1.  3  is  a 
good  example  of  this  behavior.  The  first  diffractive  minimum  (solid  line)  is 
undisturbed  moving  from  top  to  bottom.  A  maximum  (dashed  line)  and  a 
minimum  start  at  about  a  =  2.  7  and  u'  =  1.  4  (all  extrema  start  as  pairs);  they 
follow  the  diffractive  system  rather  well.  The  next  pair  of  extrema  starts 
at  about  0=  3.  9  and  u'  =  1.  6.  This  pair  is  primarily  caused  by  reflection, 
and  it  moves  more  or  less  horizontally  from  left  to  right  hand,  some  ripples 
being  superimposed.  The  pair  starting  at  about  a=  4.  7  and  u'  =  4.  1  is  part 
of  the  diffractive  system,  whereas  the  next  pair  at  a  =  5.  7  and  u'  =  4.  7  belongs 
to  a  reflective  system.  The  position  of  the  trenches  (on  the  minima)  and  the 
peaks  (on  the  maxima)  are  shown  by  large  dots.  The  dots  can  be  connected  to  ' 
distinguish  the  two  systems.  Dotted  lines  represent  the  minima  of  the  dif¬ 
fractive  system,  thin  solid  lines  the  minima  of  the  reflective  system,  and  thin 
dashed  lines  the  maxima  of  the  reflective  system. 

It  is  quite  obvious  from  this  and  other  diagrams  in  figure  3  that  trenches  appear 
always  at  the  intersection  of  the  minima  of  the  two  basic  systems  and  peaks  at 
the  intersection  of  the  maxima  of  the  two  basic  systems.  If  a  minimum  of  one 
system  intersects  with  a  maximum  of  the  other  system,  the  minima  or  the 
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Figure  2  (Concl'd) 
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Figure  4  (Concl'd) 
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"This  is  followed  by  a  careful  investigation  of  practical  limits  for  the  Mie 
region,  which  can  be  defined  only  arbitrarily.  The  lower  limit  against  the 
Rayleigh  region  is  defined  by  a  "crossover"  between  the  Rayleigh  scattering 
coefficient  and  the  Mie  scattering  coefficient.  The  upper  limit  against  the 
geometrical  optics  region  is  defined  by  the  lowest  a  value  for  which  the 
ratio  p^j  -  Pr/pm  reaches  a  value  of  less  than  25  percent.  Consequently 
in  the  geometrical  optics  region,  this  ratio  always  remains  below  25  per¬ 
cent,  For  n  =  1.33  for  example,  the  true  Mie  region  extends  from  o (  -  4.7 
to  a  =  25,  and  for  this  region,  there  is  no  other  choice  except  that  of  using 
Mie's  formulas  whereas  simple  formulas  can  be  applied  successfully  below 
and  above  these  limits. 

"Furthermore,  methods  are  developed  for  computing  the  scattering  functions 
by  simple  means.  Those  given  for  the  extended  Rayleigh  region,  and  the 
geometrical  optics  region,  lead  to  acceptable  values  for  practical  problems, 
and  errors  can  easily  be  reduced  to  less  than  10  percent.  A  new  approxi¬ 
mation  method,  derived  from  the  cross  section  theorem,  is  developed.  Since 
total  scattering  coefficients  K  are  known  to  a  very  large  extent,  this  new 
method,  which  uses  only  K  to  determine  the  scattering  functions  in  the  for¬ 
ward  direction,  can  become  of  great  value  for  tt  >  5.  Its  accuracy  is  care¬ 
fully  estimated,  and  it  is  found  that  the  errors  will  stay  below  10  percent 
with  a  few  exceptions  for  small  particle  sizes  and  will  decrease  rapidly 
with  increasing  size  of  the  aerosol. 

"Finally,  applications  of  these  theoretical  data  are  given  to  compute  the 
angular  Mie  cross  section  and  the  volume  scattering  coefficient  by  a 
graphical  method.  These  scattering  values  can  be  derived  also  for  given 
size  distributions  of  the  aerosol. " 

Corrections 

page  2,  paragraph  3,  line  2,  <  3  x  10'^cm  instead  of  <  3  x  lO^cm. 

page  12,  Ordinate,  angular  scattering  coefficient  values  for  n  =  1.  75 
revised;  see  figure,  . .  below. 

page  21,  Table  1,  extended  in  Table  I  of  Scientific  Report  7. 

page  31,  Table  II,  extended  and  revised  in  Table  II  of  Scientific  Report  7. 

page  34,  Table  III,  extended  and  revised  in  Table  III  of  Scientific  Re¬ 
port  7. 

page  38,  Table  IV,  extended  and  revised  in  Table  IV  of  Scientific  Re¬ 
port  7. 
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Reference 

n 

Reference 
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leaving  out  higher  terms  in  equation  (8)  because  terms  with  a ®  have  already 
been  left  out  in  equation  (5). 

Again  the  error  against  the  exact  solution  has  been  determined  (figures  7  and 
11  in  Scientific  Report  3),  showing  that  these  formulas  are  useful  for  a<  1  if 
n  <  2  and  X<  1. 

These  formulas,  equations  (5)  to  (9),  are  very  useful  because  this  small  size 
range  is  difficult  to  compute  using  the  exact  Mie  equations  since  the  two  recur¬ 
sion  relations  which  provide  the  input  for  the  first  term  are  losing  all  eight 
significant  figures  for  K<a)  for  large  n,  and  x  in  the  normal  programming  pro¬ 
cedure.  It  is  therefore  helpful  to  use  these  equations  for  n  and  very  small 
values  of  a,  where  special  methods  have  to  be  employed  to  obtain  reliable 
results  using  the  exact  Mie  equations. 

2.  2  FORWARD  SCATTERING 

The  forward  scattering  has  been  investigated  because  of  its  fundamental 
importance,  and  results  are  given  for  the  range  of  n  =  1.1  to  n  =  2.  Again, 
the  major  oscillations  appear  as  in  K,  and  also  the  ripples.  The  advantages 
of  various  scattering  parameters  are  outlined,  showing  the  importance  of  the 
phase  function  p^j.  The  interpolation  for  small  and  large  size  parameters  is 
possible  because  approximation  formulas  have  been  developed  for  the  occurrence 
of  the  first  "plateau,  "  its  width,  the  position  of  succeeding  plateaus,  as  well 
as  asymptotic  values  for  large  a. 

The  practical  limits  of  the  Rayleigh,  Mie,  and  diffractive  regions  have  been 
determined  as  function  of  a  and  n. 
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Figure  3  ALTITUDE  CHARTS  FOR  INTENSITY  FUNCTION  ij 

Position  of  minima  is  indicated  by  solid  lines,  of 
maxima  by  dashed  lines,  and  center  position  of 
trenches  and  peaks  by  full  circles.  Dotted  lines 
indicate  interpolation  of  diffractive  minima,  thin 
solid  lines  of  reflective  minima,  and  thin  dashed 
lines  of  reflective  maxima. 
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Scientific  Report  7,  Mie  Scattering  in  the  Forward  Area>  by  R.  Penndorf, 
AFCRL-62-1026,  RAD-TR-62-55  (September  1962),  18p. ,  reprinted  from  In¬ 
frared  Fhya.  2,  85-102  (1962). 

"The  scattering  of  light  in  the  forward  area  is  investigated  from  a  theoreti¬ 
cal  point  of  view,  and  it  is  based  on  all  available  data  computed  according 
to  the  Mie  formulas.  A  wide  range  of  refractive  indices  n  from  n  =  1.05 
to  n  =  2.0  has  been  selected.  The  behavior  of  these  scattering  coefficients 
is  investigated  and  described  in  detail.  It  is  found  that,  for  example,  while 
id  shows  major  oscillations,  only  plateaus  are  found  for  pj^j,  and  their  width 
is  related  to  the  refractive  index.  The  amplitude  of  the  major  oscillation 
decreases  with  increasing  size  parameter.  Superimposed  on  these  oscilla¬ 
tions  are  ripples,  the  phase  of  which  agrees  exactly  with  the  phase  of  the 
ripples  for  the  total  scattering  coefficient  K.  Several  ripple  systems  are 
noted  and  some  regularity  in  their  appearance  is  clearly  indicated.  The 
multi  value  ness  of  the  scattering  function  and  its  consequences  are  also 
pointed  out.  This  is  followed  by  a  careful  investigation  of  practical  limits 
for  the  Rayleigh  and  the  Mie  region. " 

Scientific  Report  8,  Approximation  Formula  for  Forward  Scattering,  byR.  Penndorf, 
AFCRL-62-1027,  RAD-TR-62-56  (September  1962),  4p. ,  reprinted  from  J.  Opt. 
Soc.  Am.  5£.  797-800  (1962). 

"By  using  the  cross  section  theorem,  an  approximation  formula  is  derived 
for  the  scattering  of  light  in  the  forward  direction.  It  is  valid  for  spherical 
aerosols  if  the  radius  of  the  particles  is  larger  than  the  wavelength  of  the 
incident  light." 

Scientific  Report  9.  Scattering  and  Extinction  Coefficients  for  Small  Absorbing 
and  Nonabsorbing  Aerosols,  byR.  Penndorf,  AFCRL-62-1 131 ,  RAD-TR-63-4 
(November  1962),  9p. ,  reprinted  from  J,  Opt.  Soc.  Am.  52,  896-904  (1962). 

"Approximation  formulas  are  derived  for  small  spherical  aerosols  (r  <X) 
by  using  the  series  expansion  of  the  exact  solution  (Mie  theory).  The  first 
three  terms  of  the  series  are  given  for  real,  complex,  and  infinite  indices 
of  refraction.  The  new  formulas  allow  for  computing  the  scattering  and  ab¬ 
sorption  coefficients  for  small  aerosols  more  simply  and  faster  than  the 
Mie  formulas.  The  error  of  our  approximation  formulas  is  determined, 
and  it  is  found  that  the  useful  size  range  is  extended  by  about  a  factor  2." 

Scientific  Report  10,  Atlas  of  Scattering  Diagrams  for  n  =  1.5,  by  R.  Penndorf, 
AFCRL-62-1 131,  RAD-TR-63-9  (January  1963),  46p. 

"Angular  scattering  diagrams  are  constructed  for  spherical  aerosols  of 
refractive  index  n  =  1.5.  They  show  the  intensity  functions  ij,  i2,  and  ij 
4  i2  as  functions  of  the  scattering  angles  for  size  parametersCt  =0.5(0.5)10. 


-25- 


I.  INTRODUCTION 


This  bibliography  is  based  on  notes  collected  over  the  past  10  years,  and  it 
summarizes  what  is  known  about  existing  numerical  computations  using  the  Mie 
theory.  The  computations  for  single  scattering  are  only  considered,  and  those 
for  multiple  scattering  have  been  excluded. 

Lists  already  exist;  for  example,  some  are  given  by  van  de  Hulst  in  his  book  on 
pages  165-171  and  273-275,  and  another  by  Penndorf  (reference  42  in  section 
3.  3  of  this  appendix).  However,  these  and  similar  lists  and  bibliographies  are 
superseded  by  the  vast  amount  of  new  computations,  which  are  a  consequence  of 
the  electronic  computers.  Hence,  there  seems  to  be  a  need  for  an  updated  list. 

The  user  of  this  bibliography  should  be  familiar  with  the  field.  That  means  he 
should  have  looked  in  van  de  Hulst' s  book,  and  he  should  be  familiar  with  the 
basic  equations;  otherwise,  the  tables  will  be  useless  to  him.  Unfortunately, 
no  general  agreement  exists  on  symbols,  and  this  makes  it  difficult  for  the  user 
of  tables  scattered  over  a  large  number  of  journals.  He  should  in  all  cases 
read  very  carefully  the  definitions  of  each  author,  and  he  should  compare  them 
with  those  used  in  van  de  Hulst's  book  to  derive  the  proper  conversion  factors. 

It  seems  hopeless  to  indicate  all  the  differences  and  conversion  factors  in  the 
following  tables. 

For  reasons  of  convenience,  reference  alwayB  iB  made  to  van  de  Hulst's  book. 

It  is  an  excellent  book,  and  it  is  available  in  all  libraries;  whereas,  less  known 
references  may  be  harder  to  obtain. 
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II.  SCIENTIFIC  RESULTS 


In  this  section,  a  short  description  of  the  highlights  of  the  results  will  be 
attempted.  Three  areas  have  been  investigated:  (a)  the  total  scattering  coef¬ 
ficients,  (b)  the  forward  scattering,  and  (c)  the  angular  scattering  coefficients. 

2.  1  TOTAL  SCATTERING  COEFFICIENT 

The  total  scattering  coefficient  K  has  been  investigated  to  find  an  approximation 
method  which  allows  to  construct  a  curve  of  K  versus  a,  the  size  parameter, 
for  any  refractive  index  n  <  2.  Maxima  and  minima  of  K  occur  at  specific  size 
parameters.  Choosing  a  normalized  size  parameter  p-  2  a(n  -  1),  it  was 
found  that 


py  =  2ff(y  +  1/4)  (1) 

is  a  good  approximation;  the  subscript  y  stands  for  the  order  of  the  extreme 
values,  y  =  1 ,  2,  ...  ;  the  +  sign  applies  to  the  maxima,  the  -  sign  to  the 
minima.  For  n  =  1,  equation  (1)  is  correct;  but  for  n  >  1,  a  correction  term 
should  be  added,  so  that 

py  n  =  2ff(y+  1/4)+  0.  3  (n-  1)  (2) 

represents  the  best  approximation. 


The  absolute  value  of  Ky  at  the  extreme  values  has  also  been  obtained.  lt_was 
found  that  a  resonably  good  approximation  for  a  smoothed  extreme  value  Ky 
can  be  described  by 


4  29M 

'V2*  Py 


51 M 


(3) 


for  the  maxima,  and 


-  4  4  8.  01 M  27.  3M 

K"  =  2 - -  - - -  (4) 

py  py2  py  py2 


for  the  minima.  Here  M  =  (n^  -  1 )  /  (n^  +  2).  A  graphical  method  to  determine 
smoothed  Kfor  any  aand  any  n  has  been  described,  based  on  computations  of 

py,  K'y,  and  K'^  only. 


t 
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page  41 ,  Figure  16,  the  area  1 . 6  <  n  <  2.0  can  be  revised  based  on 
new  data. 

page  42,  Figure  17,  the  curve  for  n  =  1.5  can  be  revised  for  n  >  30 
using  the  publications  by  Giese,  Bullrich,  and  de  Barry  (1962). 

page  53,  Figure  19,  arrow  on  right-hand  side  goes  in  wrong  direction. 

Scientific  Report  3,  Scattering  Coefficients  for  Absorbing  and  Nonabsorbing 
Aerosols,  by  R.  Penndorf,  AFCRL-TN-60-667,  RAD-TR-60-27  (October  1960), 
92p. 


"Applied  problems  involving  light  scattering  by  aerosols  demand  reliable 
numerical  data  of  the  scattering  and  absorption  coefficients.  When  the  size 
of  the  aerosol  is  smaller  than  the  wavelength  of  the  incident  light,  approxi¬ 
mation  formulas  can  be  derived  from  the  Mie  theory.  The  asymptotic  ex¬ 
pansion  of  the  exact  solution  is  given,  where  the  leading  term  is  identical 
with  the  Rayleigh  approximation.  For  the  first  time,  the  first  three  terms 
as  function  of  the  size  parameter  a  are  correctly  derived  for  nonabsorbing, 
absorbing,  and  metallic  spherical  particles.  Next,  the  error  of  our  approxi¬ 
mation  is  determined  using  one,  two,  and  three  terms.  This  approximation 
extends  the  useful  size  range  by  about  a  factor  2  for  nonabsorbing  aerosols. 

"The  new  approximation  allows  for  computing  the  scattering  and  absorption 
coefficients  for  small  aerosols  simpler  and  faster  than  the  Mie  formula. 
Numerical  results  are  given  for  nonabsorbing  aerosols  having  refractive 
indices  between  n  =  1  and  n  =  2,  For  absorbing  spheres,  the  numerical 
values  are  given  for  the  coefficients  in  the  appropriate  formulas.  A  wide 
range  of  indices  of  refraction  n  =  n  -  ix  are  selected,  values  between  1  and 
9  are  chosen  for  n  and  between  1  and  4  for  x. 

"A  comparison  is  made  between  results  obtained  by  the  approximation 
formulas  and  the  Mie  formulas  for  ff  =  1.25  -  ix,  7T  =  1.29  -  ix,  ff  =  1. 50  -  ix, 
and  IT  =  1.75  -  ix.  The  new  approximation  formulas  lead  to  errors  in  the 
total  extinction  coefficient  K(e)  of  less  than  20  percent  if  a  <0.8  and  x  <  1. 
The  ratio  K(*)/K(e)  shows  that  our  formulas  can  be  used  up  to  a  =  1.0  be¬ 
cause  the  error  resulting  from  the  approximations  remains  small  and  within 
tolerable  limits. 

"Finally,  numerical  data  are  computed  for  IT  =  1.29  -  ix,  IT  =  1.33  -  ix,  and 
n  =  1.40  -  ix.  For  the  case  ?T  =  1.29  -  ix,  a  complete  new  set  of  K(e)  values 
has  been  derived  to  a  =  7  by  carefully  interpolating  those  new  data  as  well 
as  other  existing  data." 
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8.  Footnote! .  At  times  it  seemed  necessary  to  add  a  footnote  to 


a  specific  listing.  This  is  done  only  rarely  and  to  save  space  in  the  last 
coluan.  The  footnotes  are  collected  at  the  end  of  each  section,  (3*1,  ^*1> 
5.1  and  6.1). 

Footnotes  are  indicated  by  (FI,  12...)  in  the  Inst  column. 
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other  than  In  connection  with  a  definitely  related 
government  procurement  operation,  the  U.  S. 
Government  thereby  Incurs  no  responsibility,  nor  any 
obligation  whatsoever;  and  the  fact  that  the  Govern¬ 
ment  may  have  formulated,  furnished,  or  In  any  way 
supplied  the  said  drawings,  specifications,  or  other 
data  Is  not  to  be  regarded  by  Implication  or  other¬ 
wise  as  In  any  manner  licensing  the  holder  or  any 
other  person  or  corporation,  or  conveying  any  rights 
or  permission  to  manufacture,  use  or  sell  any 
patented  Invention  that  may  In  any  way  be  related 
thereto. 
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Figure  5  ALTITUDE  CHARTS  FOR  INTENSITY  FUNCTION  i2 

Position  of  minima  is  indicated  by  solid  lines,  of 
maxima  by  dashed  lines,  and  center  position  of 
trenches  and  peaks  by  full  circles.  Dotted  lines 
indicate  interpolation  of  diffractive  minima,  thin 
solid  lines  of  reflective  minima,  and  thin  dashed 
lines  of  reflective  maxima. 
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VI.  APPENDIXES 

6.  1  BIBLIOGRAPHY  OF  NUMERICAL  COMPUTATIONS  ON 
SCATTERING  AND  ABSORPTION  OF 
ELECTROMAGNETIC  RADIATION 
FOR  SPHERICAL  PARTICLES 
BASED  ON  THE  MIE  THEORY 

6.  2  ATLAS  OF  SCATTERING  DIAGRAMS  FOR  n  =  1.  1,  1.2, 
1.  33,  1.4,  1.  44,  AND  1.  5 
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TABUS  IV 


* 


ABBREVIATION 

Tb*  following  abbreviation*  are  used  in  all  tables  listing  the  existing 
computations  (sections  3*1,  ^.1,  5*1  and  6.1).  The  abbreviations  are  se¬ 
lected  on  the  basis  of  avoiding  subscripts,  to  facilitate  typing  and  to  save 

space. 
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III.  NONABBCRBIWO  SPHERES 


3.1  List  of  existing  computations  for  affl,  bm,  S1,  Sg,  1^  lg  and  K  for  mono- 

disperse  isotropic  spherical  particles  (Dielectrl.c  and  perfectly  con¬ 
ducting  particle*). 

The  following  table  requires  some  explanation.  The  first  column  lists 
a  running  number.  This  number  agrees  with  the  number  in  the  references, 
section  3*3*  The  second  column  lists  the  "name  of  the  author (s)  and  the 
year  of  publication" .  The  symbol  U  after  the  year  means  that  the  data  are 
unpublished.  The  third  column  contains  the  refractive  index  n. 

The  fourth  column  lists  the  "range  of  values" .  The  accepted  system 
is  used;  1(1)10(10)  100  means  steps  of  /\CX  *=  1  cure  used  between  <X*  1  and 
<X  «  10  and  steps  of  A0<  m  10  axe  given  between  10  and  100.  In  the 
older  literature  we  find  that  uneven  steps  of  (X  have  been  taken  and  therefore 
we  report  the  results  in  the  following  system;  O.5-5.7  £l8j  which  means 
the  lowest  CK value  is  0.5,  the  highest  is  5.7  and  18  0( values  between  these 
limits  are  listed. 

The  fifth  column  lists  the  "tabulated  quantities".  This  is  the  most 
difficult  column  to  devise  and  a  simplified  system  similar  to  that  of  v.d. 
Hulst  and  Fenadorf  is  used.  To  use  single  space  typing  all  subscript  are 
avoided  If  possible.  The  order  of  symbols  (see  here  especially  Table  IV 
for  abbreviations  is  a,  b,  then  S,  then  K,  then  1  or  I,  followed  by  less 
frequent  symbols,  such  as  R,  P.  The  listing  "i  0  »  0(10)(l80)"  means  i^  and 
12  are  listed  for  scattering  angle  0  -  0°  to  180°  in  steps  of  A®  ■  10°. 
Without  exception,  0  ■  0  represents  forward  scattering.  Irrespective  of  what 
the  author  has  used  as  definition  for  the  forward  angle.  This  is  done  so 
that  the  reader  looking  for  a  specific  scattering  angle  can  use  the  table, 
fbs  syAol  •  in  the  next  column  has  to  be  watched.  It  Indicates  that  the 


Footnotes 


#1  Information  copied  from  v.d.  Bulat.  Seems  doubtful.  I  believe  Rayleigh 
used  radius  y*  not  P<  and  log  cos  0.  Mo  thorough  check  has  been  made  be¬ 
cause  the  results  have  only  historic  Interest. 

#2  It  is  not  quite  clear  vbat  refractive  index  has  been  used.  In  the  first 
paper  n"  =  ^  (1-0.031)  is  mentioned,  in  the  other  papers  no  value  for  n 
is  given.  I  believe,  however,  that  a^  ag,  a^,  b1,  b£  are  Identical  In 
all  three  papers.  The  value  of  CX Is  not  correctly  stated  In  v.d.  Bulst 
(p.  168 ),  because  Ruedy  used  CX  ’  ■  2  77^  . 
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V.  HETERODISPEBSE  SY8TEMB  OF  ISOTROPIC  SPHERICAL  PARTICLES 


5.1  List  of  existing  computations 

The  list  has  been  prepared  in  a  form  similar  to  3>1  and  4.1.  The  important 
column  here  is  the  "size  distribution"  used  by  the  author.  We  have  only  indi¬ 
cated  the  type  of  distribution  function  used;  for  detail  one  has  to  read  the 
original  paper.  It  seems  impossible  to  list  an  exact  formula  for  each  case. 
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6.2I+  Atlas  of  scattering  diagrams 
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SCATTERING  ANGLE  6,  IN  DEGREES 


Figure  2  DIAGRAMS  OF  ANGULAR  POSITION  OF  MAXIMA  (BRIGHT  RINGS) 
AND  MINIMA  (DARK  RINGS)  FOR  INTENSITY  FUNCTION  ij 

Po«ition  of  minima  is  indicated  by  solid  lines,  of  maxima  by 
dashed  lines,  and  center  position  of  trenches  and  peaks  by 
full  circles.  Dotted  lines  indicate  interpolation  of  diffractive 
minima,  thin  solid  lines  of  reflective  minima,  and  thin  dashed 
lines  of  reflective  maxima 
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"The  basic  data  have  been  computed  using  Mie's  theory,  and  a  graphical  in¬ 
terpolation  technique  (altitude  chart  technique)  has  been  used  to  determine 
accurately  the  position  and  numerical  value  of  the  maxima  and  minima. 

"The  scattering  diagrams  presented  in  this  atlas  allow  to  determine  the 
intensity  functions  ij ,  i£>  and  ij  i 2  for  any  desired  scattering  angle  with 
a  high  degree  of  reliability.  The  basic  data  limit  the  construction  of  re¬ 
liable  scattering  diagrams  to  a  <  10." 

Corrections 

Figures  3-1  and  3-4  have  been  inadvertently  interchanged, 
page  9,  Figure  3-4,  instead  of  Figure  3. 1. 
page  13,  Figure  3-1,  instead  of  Figure  3-4. 
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4.2  Li»t  of  Coigutatioaa  arrsaged  according  to  Refractive  Index 

Instead  of  a  listing,  we  have  prepared  a  figure  la  which  the  ordlaate 
represents^,  the  real  part  aad  the  abscissa  *,  the  lwegi nary  part.  Each 
point  la  the  plaae  represents  then  a  refractive  Index  a'  for  which  the  sane 
type  of  computation  has  been  accomplished. 
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2.  3  ANGULAR  SCATTERING  COEFFICIENTS 


A  large  part  of  the  effort  has  been  devoted  to  an  interpretation  of  the  angular 
scattering  coefficients  ij  and  for  the  refractive  indices  n  =  1.  33  to  1.  5.  The 
objective  has  been  to  compose  an  atlas  of  angular  scattering  coefficients, 
which  is  now  completed  and  given  as  appendix  6.  2.  Since  the  basic  data  are 
computed  in  steps  of  5  degrees  (for  n  *  1.33)andl0  degrees  (for  all  other  refrac¬ 
tive  indices),  interpolative  methods  had  to  be  devised  so  that  continuous  curves 
for  the  scattering  coefficients  could  be  drawn.  It  was  found  that  the  "altitude 
chart"  offers  the  best  method  for  graphically  interpolating  scattering  data.  It 
is  important  to  find  systems  for  the  maxima  (bright  rings)  and  minima  (dark 
rings)  because  the  interpolation  of  the  computed  data  is  very  much  simplified 
if  the  angular  position  and  the  absolute  value  of  the  extreme  values  of  ij  and 
i2  are  know. 

Many  different  and  laborious  attempts  have  been  made  to  find  a  systematic 
behavior  of  the  maxima  and  minima  and  to  understand  it.  All  the  data  have 
been  investigated  again,  and  the  present  position,  modifying  earlier  results 
reported  in  Scientific  Reports  4,  5,  and  10,  follows.  It  is  believed  now  that 
the  first  interpretation  of  the  systematic  behavior  of  the  bright  and  dark  rings 
is  not  always  correct. 

For  graphical  interpolation,  the  "altitude  chart"  method  has  been  developed 
for  ij,  i£»  and  ig  and  described  in  Scientific  Reports  5  and  10.  Based  on  the 
recent  experience  in  finding  a  system  for  the  diffractive  as  well  as  the  reflec¬ 
tive  minima,  it  would  seem  that  this  altitude  chart  can  be  improved.  The  new 
altitude  chart  is  based  on  the  two  parameters  u  =  a  Bin  0  and  a,  instead  of  u 
and  p  .  Hence,  one  can  use  the  symmetry  at  8  =  90  degrees  and  have  two 
separate  areas;  i.  e. ,  one  for  0  =  0  to  90  degress,  and  the  other  for  0  =  90  to 
180  degrees;  at  0  =  90  degrees  the  two  systems  become  identical.  Ordinate 
a  and  abscissa  u  =  a  sin  0  are  used,  covering  all  angles  between  0=0  and 
90  degrees.  At  0  =  90  degrees,  the  system  is  changed  so  that,  for  the  realm 
0  =  90  to  180  degrees,  the  abscissa  is  aand  the  ordinate  u'  =  or  sin  0',  where 
0'  =  ff-  0L  In  figure  1,  a  schema  of  the  new  altitude  chart  is  shown.  It  will 
be  seen,  in  the  examples  which  follow,  that  the  diffractive  minima  and  maxima 
are  starting  now  at  the  top  and  move  downward,  whereas  in  the  old  system  they 
start  at  the  right-hand  side,  move  horizontally  to  the  left-hand  side,  a  knee 
appears  at  0  =  90  degrees,  and  only  in  the  forward  scattering  area  do  they  move 
downward.  The  reflective  minima  will  now  move  more  or  less  horizontally 
from  left  to  right  hand,  whereas  in  the  old  system  each  one  possesses  a 
different  inclination  to  the  basic  parameters  pand  u.  It  is  believed  that  this 
new  altitude  chart  is  better  than  the  old  chart  in  recognizing  the  two  basic 
systems,  especially  around  0  =  90  degrees,  in  aiding  the  interpolation  and  in 
determining  the  position  and  depth  of  the  minima  at  points  of  interference  of 
the  diffractive  and  reflective  systems. 
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The  difference  between  succeeding  minima  is  then  given  by 

Act  =—  (IT  -  1)  -  I  +  1  (n  -  l)2 
n  4  4  4 

=  ;(ff-;)-;[1'(n'1)2]*  (12) 

These  approximation  formulas  are  valid  only  for  n  <  1.5;  they  have  not  been 
checked  for  n  >  1.5. 

In  summary,  graphical  methods  for  interpolation  of  computed  data  have  been 
developed  which  promise  to  lead  to  reliable  interpolation  for  practical  applica¬ 
tion  of  the  Mie  theory. 
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Scientific  Report  5,  Atlas  of  Scattering  Diagrams  for  n  =  1.33,by  R.  Penndorf, 
AFC RL- 1044,  RAD-TR-61-32  (October  1961),  69p. 

"Angular  scattering  diagrams  are  constructed  for  spherical  aerosols  of  re¬ 
fractive  index  n  =  1.33.  They  show  the  intensity  function  ij  and  i 2  as  func¬ 
tion  of  the  scattering  angle  for  size  parameters  a  =  0.5  (0.5)  15.  The  basic 
data  have  been  computed  using  Mie's  theory  and  graphical  interpolation 
techniques  have  been  used  to  fill  the  gaps. 

"The  diagrams  presented  in  this  atlas  allow  to  determine  the  intensity  func¬ 
tions  ij  and  i2  for  any  desired  scattering  angle  with  a  high  degree  of  re¬ 
liability.  " 

Corrections 

Figure  1,  the  insert  is  taken  from  Giese  for  n  =  1.33  and  a  =  10. 

Figure  2,  a  new  schema  of  an  improved  altitude  chart  is  discussed  in 
section  II  of  this  final  report. 

Figure  9,  a  new  interpretation  of  this  figure  can  be  given  based  on  our 
discussion  in  section  II  of  this  final  report. 

page  24,  line  3,  polynomials. 

page  25,  Table  I,  the  "number  of  rings"  given  under  "Penndorf"  is  now 
considered  too  large.  Since  some  of  the  minima  can  "vanish, "  it  is 
believed  that  the  number  of  rings  for  a  >  12  is  overestimated.  The 
"suggested  steps  of  A0  in  degrees,"  however,  will  not  change  very 
much  because  the  interval  between  minima  can  still  be  on  the  low  side 
as  indicated. 

page  29,  it  should  be  added  that,  in  the  atlas,  ij  is  shown  by  a  solid 
line  and  i2  by  a  dashed  line. 

page  56,  Dr.  Diermendjian  (personal  communication)  informed  us  that 
his  values  for  n  =  1.33,  0  =  45  degrees,  and  rt  =  13,  14,  and  15  differ 
somewhat  from  those  given  in  the  report.  The  differences  have  not 
been  settled  as  yet. 

Scientific  Report  6,  Angular  Mie  Scattering,  by  R.  Penndorf,  AFCRL-62-1025, 
RAD-TR-62-54  (September  1962),  7p. ,  reprinted  from  J.  Opt.  Soc.  Am.  5£, 
402-408  (1962). 

"The  basic  Mie  formulas  are  explained,  and  several  useful  scattering  func¬ 
tions  are  defined,  such  as  the  angular  Mie  scattering  coefficients  ig  and  the 
Mie  phase  function  pM.  Finally,  the  scattering  pattern  for  the  forward,  side¬ 
ward,  and  backward  area  is  surveyed,  based  on  data  for  the  refractive  in¬ 
dex  n  =  1.33." 
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plex.  The  formulas  have  been  derived  from  Mie's  formula  using  the  series 
expansion  method.  Their  limits  have  been  described." 

Penndorf,  R. ,  Approximation  formula  for  forward  scattering,  J.  Opt.  Soc.  Am. 
52,  797-800  (1962). 

See  under  Scientific  Report  8  for  abstract. 

Penndorf,  R. ,  Mie  scattering  in  the  forward  area.  Infrared  Fliys.  2,  85-102 
(1962). 

See  under  Scientific  Report  2  for  abstract. 

Penndorf,  R. ,  Scattering  and  extinction  coefficients  for  small  absorbing  and 
nonabsorbing  aerosols,  J.  Opt.  Soc.  Am.  52,  896-904  (1962). 

See  under  Scientific  Report  9  for  abstract. 

Penndorf,  R. ,  Scattering  Diagrams  in  the  Mie  Region,  Conf.  Electromag. 

Scat.  (ICES),  p.  73-86  (1963). 

"Scattering  diagrams  showing  the  intensity  function  ij  as  function  of  the 
scattering  angle  have  been  constructed  for  4  refractive  indices  between 
n  =  1.  33  and  1.50.  They  are  given  in  an  atlas  for  size  parameters  a  =  0.  5 
(0.  5)  10  and  in  the  case  of  n  =  1.  33  to  a-  15.  The  basic  data  have  been  conv- 
puted  using  Mie's  theory,  and  graphical  interpolation  techniques  have  been 
used  to  fill  the  gaps.  The  various  interpolation  techniques  are  discussed  in 
detail,  especially  the  altitude  chart  techniques.  They  are  very  important  to 
fill  the  gaps  in  the  basic  data. 

"The  scattering  diagrams  allow  determination  of  the  intensity  functions  ij 
and  for  any  desired  scattering  angle  with  a  high  degree  of  reliability. 
Examples  of  the  results  are  shown. " 
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so  entrenched  In  the  radar  literature  that  it  Is  hopeless  to  change  It. 
Equation  3.3  (Table  III)  and  van  de  Hulst  (p.284)  show  the  relationship 
between  cr  and  1^  and  ig.  One  has  to  be  careful  if  the  author  uses  polar¬ 
ized  or  unpolarized  radiation  for  the  source  or  receiver,  or  both.  Hence, 
a  factor  2  can  appear.  Again  the  definitions  of  the  author  have  to  be 
watched  before  using  tabulated  data. 

6.  Letters  and  Symbols.  No  standard  notation  exists.  Hence, 
the  reader  will  find  all  kinds  of  letters  for  the  sauna  quantity.  It  scene 
hopeless  to  include  such  a  list,  although  it  can  be  prepared. 

for  exanple,  for  the  amplitude  functions  we  find  .2^, 

#  # 

ij,  lgj  E p,  E e;  and  Aj,  Ag.  A  complete  list  of  conversion  factors, 
although  extremely  helpful,  requires  a  lot  of  work  but  nay  be  needed  only 
by  a  few. 

Again  we  can  only  warn  the  reader  to  pay  attention  to  the  defini¬ 
tions  used  by  each  author  and  compare  then  with  Table  I,  II,  and  III  or 
van  da  Bulat's  book.  FTon  there  the  reader  will  be  able  to  arrive  at  a 
conversion  factor,  if  it  is  really  necessary. 

7.  Odd  parameters.  Problems  arose  in  compiling  our  tables  when 
the  author  did  not  use  our  basic  system  of  n  and  c*  in  giving  his  results, 
for  example,  we  will  find  that  an  author  specified  only  A  and  T"  and  refers 
to  a  specific  cheaioal  substance  for  the  refractive  index  at  those  A 'a. 
frequently  the  eouroe  for  n  ie  listed,  but  not  the  adapted  value  of  n. 

We  have  not  checked  further  since  we  are  not  sure  what  values  have  been 
used  for  the  cogitations.  In  such  cases  we  give  the  information  as  it 
appears  in  the  listed  source  and  have  refrained  fron  a  conversion  into  our 
fiysten.  This  applies  especially  to  section  IV. 
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The  applications  of  our  results  are  manifold,  and  since  scattering  by  spherical 
particles  has  been  proven  as  an  extremely  helpful  tool  in  many  disciplines,  no 
listing  of  possible  applications  of  the  results  will  be  given. 
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Figure  1  SCHEMA  OF  MODIFIED  ALTITUDE  CHART 

Scattering  angles  (9  -  0  to  180  degrees)  are 
represented  by  straight  lines.  Data  are 
plotted  at  intersection  of  a  with  9  . 
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Figure  4  DIAGRAMS  OF  ANGULAR  POSITION  OF  MAXIMA  (BRIGHT  RINGS) 
AND  MINIMA  (DAT.K  RINGS)  FOR  INTENSITY  FUNCTION  i2 

Position  of  minima  is  indicated  by  solid  lines,  of  maxima  by 
dashed  lines,  and  center  position  of  trenches  and  peaks  by 
full  circles.  Dotted  lines  indicate  interpolation  of  diffractive 
minima,  thin  solid  lines  of  reflective  minima,  and  thin  dashed 
lines  of  reflective  maxima. 
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maxima  are  shallow.  This  is  also  the  position  where  a  system  may  suddenly 
"vanish,  "  or  where  an  analysis  of  the  minima  suddenly  jumps  from  one  system 
into  the  other  system.  Theoretically,  this  can  be  explained  on  the  basis  of 
ray  optics.  At  the  intersection  of  the  two  systems,  both  rays  are  of  opposite 
phase,  so  that  the  intensity  becomes  extremely  low.  For  the  maxima,  the 
two  systems  must  be  in  phase  so  that  the  intensity  reaches  peaks.  It  seems 
possible  to  investigate  this  phase  effect  for  some  of  the  points  and  corroborate 
intuition.  The  cases  for  n  =  1.  33  and  n  =  1.  486  are  the  best  examples  for 
this  interpretation  although  only  the  minima  are  shown. 

For  the  intensity  function  i£i  the  conditions  are  somewhat  different  because 
the  minimum  around  0  =  90  degrees  exercises  a  dominant  influence.  This 
situation  presents  itself  already  for  r.  =  1.  05  (figures  4  and  5).  The  dif¬ 
fractive  minima  and  maxima  can  be  interpolated  around  0  =  90  degrees  (dotted 
lines),  and  they  cross  each  other  where  they  intersect  the  90-degree  min¬ 
imum.  The  90-degree  minimum  moves  always  into  the  backward  area 
where  it  intersects  the  diffractive  minima.  Some  of  this  behavior  is  seen 
better  in  figure  4  than  in  figure  5.  For  n  =  1.  15,  the  matter  starts  to  look 
more  complicated  because  the  second  and  third  diffractive  minimum  does 
not  start  at  0=180  degrees, but  a  pair  of  maxima  and  minima  starts  around 
135  and  150  degrees.  The  case  for  n  =  1.  20  (figures  4  and  5)  shows  the  con¬ 
dition  very  well.  The  minimum  starting  at  0  =  90  degrees  for  a=  0.  1  under¬ 
goes  large  angular  fluctuations  between  0  =  100  and  140  degrees;  its  behavior 
is  best  seen  in  the  altitude  chart,  where  it  follows  a  diffractive  minimum 
for  a  while  and  then  moves  suddenly  to  large  scattering  angles  until  it  meets 
the  next  diffractive  minimum.  The  first  reflective  system  is  indicated  by 
these  sharp  tips.  Inflective  points,  instead  of  well  developed  minima,  are 
shown  by  the  sign  .  The  altitude  chartfor  n  =  1 .  3  allows  to  find  the  first  and 
second  reflective  system  very  clearly,  and  figure  5  shows  the  solution  for  n=  1.  33 
and  n  =  1.486  as  the  best  examples  for  the  existence  of  these  two  systems.  The 
minima  for  i ^  and  small  a,  say,  CK  below  6,  are  very  often  interrupted,  they  just 
"disappear"  on  figure  4;  but  the  altitude  charts  reveal  clearly  that  two  systems 
exist,  and  the  position  of  the  minima  can  be  interpolated  without  difficulty. 

On  the  various  diagrams,  the  diffractive  minima  have  been  interpolated  (dotted 
lines)  through  the  backscattering  area  to  6  =  iso  degrees.  The  value  for 8  = 

180  degrees  has  been  computed  according  to  the  following  scheme.  It  was  in-  ' 
vestigated  whether  a  system  existed  for  the  beginning  of  the  diffractive  minima 
at  9  =  180  degrees.  Finally,  such  a  system  emerged,  so  that  for  ij  and  i 2 

aj=f  [n+j  <£-2n+n2>]’ 

where  tt.  denotes  the  a  value  at  which  a  new  diffractive  minimum  begins  at 
0  =  180  degrees,  and  j  =  1,2,3  stands  for  the  first,  second,  third,...,  diffrac¬ 
tive  minimum.  Thus,  one  is  able  to  draw  the  "undisturbed"  curve  for  each  dif¬ 
fractive  minimum  (dotted  line).  A  system  must  also  exist  for  the  reflective 
minima,  but  so  far  an  approximation  formula  has  not  been  found. 
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II.  BESGRIPTCOS  OF  SIMBOI£  AMD  TAMES 


To  understand  the  tables  listed  In  this  report,  the  following  quanti¬ 
ties  and  formulas  are  needed.  We  have  excluded  a  discussion  and  listing 
of  all  auxiliary  functions,  such  as  Bessel  functions  and  Legendre  poly¬ 
nomials;  some  of  these  can  be  found  on  p.  165  in  v.d.  Bulat's  book,  al¬ 
though  his  list  is  incomplete. 

The  definitions,  symbols,  and  formula  are  arranged  in  tabular  form, 
similar  to  Havksley's.  This  seems  to  be  the  most  concise  form.  Table  1 
contains  the  definitions  and  symbols.  Table  II  the  formula  for  total  ex¬ 
tinction  and  scattering.  Table  III  the  formula  for  angular  scattering,  and 
Table  IV  the  abbreviations  used  in  sections  III  to  V. 

Some  notes  and  explanations  are  deemed  advisable. 

1.  The  refractive  index  is 

IT  -  n  -  ilC  .  (1) 

This  Is  the  complex  index  of  refraction;  for  a  nonabsorbing  particle  If  ■  0. 
Some  authors  prefer  the  sysfcol  m  Instead  of  if.  So  far  no  agreement  exists 
about  the  use  of  m  or  n.  The  SOI  Comnisslon  (Physics  to-day,  1962)  suggest 
the  symbol  n.  Ve  should  also  note  that  some  authors  prefer 

n  -  n  (1  -  IK).  (2) 

In  addition,  the  dielectric  constant 

L-L'-it"  (3) 

Is  given,  especially  in  radar  work.  This  can  be  converted  into  n  by  stand¬ 
ard  procedure  (  f  s  ■**’)• 

2.  The  total  Hie  coefficients  K,  K*,  K*,  if  pose  no  problem,  ex¬ 
cept  K/2  Is  given  in  the  older  literature. 

3.  The  angular  scattering  coefficients  are  1^,  ig.  The  scatter- 
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author  has  used  another  definition  for  forward  scattering,  mostly  O'  -  100°. 
In  those  cases  where  only  i^  is  given  it  has  been  spelled  out.  If  com¬ 
pletely  different  quantities  are  listed  like  the  angular  distribution  coeffi 
dent  (references  48,  49),  the  quantity  has  been  spelled  out. 

The  sixth  coluac  lists  the  form  of  "data  presentation'’ .  First  the 
symbol  g  indicates  that  the  definition  O'  -  l8o  far  forward  scattering  is 
used  by  the  author;  hence,  first  the  tables  for  and  lg  have  to  be  watched 
and  secondly  the  tables  for  a^  and  ba  have  to  be  converted  according  to 
equation  4.  The  systole  T  and  0  tell  you  whether  the  author  gives  the  re¬ 
sults  in  tabular  or  graphical  form  or  both.  Finally,  the  letter  H  indicates 
that  the  source  has  not  been  checked.  This  applies  to  reports  which  I  have 
not  seen,  but  the  results  have  been  quoted  from  another  source,  which  I  be¬ 
lieve  to  be  reliable. 

The  seventh  and  last  column  Indicated  the  "type  of  computation", 
namely,  whether  a  desk  computer  was  used  or  an  electronic  computer  or  any 
other  type.  This  might  be  helpful  because  the  program  nay  be  available 
to  thoae  having  similar  machines  at  their  disposal.  The  older  computations 
are  naturally  all  "hand-aada"  or  by  a  simple  desk  computer  and  liable  to 
errors  of  eosgutatlon.  In  those  computations  the  tern  "reliable"  or  "un¬ 
reliable”  Indicate  what  we  know  about  the  result  in  order  to  warm  the  user. 
The  term  unreliable  comes  19  whenever  the  author  did  not  use  enough  tern  of 
the  series  (equation  2.$)*  The  electronic  counter  should  lead  to  correct 
results,  but  errors  can  creep  in,  namely  the  program  may  contain  a  alight 
error,  or  the  machine  used  an  old  tape  at  a  certain  portion  of  the  coaputa- 
tion  or  the  programs  were  written  for  one  machine  but  used  by  another 
machine.  Those  errors  are  hard  to  find  by  Just  looking  at  tables;  however, 
graphical  representations  will  show  it.  The  content  and  accuracy  of  tabu- 
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Author  (a) 

Refractive 

Data 

No. 

year  of  Publ. 

Index 

Rang*  of  o(  Value* 

Tabulated  Quaatltie* 

Preaentatlon 

■eeesjKMI 

13 

Fuedy 

1.33 

<X4  6 

K  (#),a,b 

0 

D,  F(2) 

1943/44 

370^/4,  3^8,^/2 

II  0-0-180 

0 

D,  rs 

33/4,  5F- 

Jp  0-50-130 

0 

14 

Langmuir 

1.50 

1-10  Gs4j 

1/2K 

T 

rel;  D,  4a 

1943 

2.0 

1-5  [21] 

1/2K 

aee  ref. (19) 

15 

Rubin* tain 

Pel  lam 

1943 

OO 

0-10 

R  0-180 

N 

n  ! 

16 

Peary,  Scott 
1944 

0-25 

R  0-180 

17 

van  de  Hulat 

1.0 

?<  23  rl 

K 

T 

rel;  D  | 

1946 

90 

0.5-10  15J 

K 

i1 

| 

18 

Ho  11 

4/3 

0.3-18  Oto] 

1/2K,  a,  b. 

T,G  • 

rel;D;  4,5p 

1946/48 

0.3-6  [16] 

1,F  0-0(10)180 

T,G  • 

4p 

6.6-18  [24] 

1,  0-0(90)180 

T  • 

4p 

7T,  £0-0(10)180 

T 

4p 

19 

Barnea,  Kenyon 

Zalaer,  LaMar 

1947 

1.50 

0.5-12  fee] 

K,  Hating  of  ref. 

(14,23) 

T,0 

rel;  D,  4a 

1  20 

Hyde 

1.75 

R  0-180 

1947 

21 

Spendley 

1947 

OO 

0-12 

R  «*180 

1 

22 

Houghton 

4/3 

9-24  [50] 

K 

T,0 

rel;  ex- 

C talker 

trmpol.  to 

1949 

CX-  50;  D 

23 

Lowan 

1.33 

0. 5-6.0  [15] 

a,b,K,  1  0-0(10)180 

•.  r 

rel;  D, 

1949 

1.44 

Mat 

3*(K) 

1.50 

MM 

4*  or  6p 

2.0 

6.5-12  [lo] 

(*,b,l) 

1.50 

K 

1.44(0.1) 

1-55 

0. 5-7.0  po] 

K 

24 

Riley 

1949 

1.486 

0.5(0. 1)3 

a,b,K,i  0-0(5)180 

9 

rel;  D 

48- 


73 


VI.  TWO  OR  MORE  CONCENTKIC  SPHERES  OF  DIFFERENT  REFRACTIVE  INDEX 


6.1  List  of  computations 

The  11 8 t  has  been  prepared  in  a  form  similar  to  3*1  and  4.1.  The  re¬ 
fractive  indices  n.  and  n  indicate  the  refractive  index  for  the  inner  core 
i  a 

and  the  outer  coating,  similarly  the  symbols  (X.  and  o(  or  r.  and  r 

l  a  l  & 

indicate  the  size  parameter  or  radius  for  the  core  and  the  coating; 
i.e.,  ra  is  the  thickness  of  the  center  shell;  q  =  r^/ (r^  +  ra). 
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INTENSITY  FUNCTION 


102 
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INTENSITY  FUNCTION 


9 


131 


INTENSITY  FUNCTION 


142 


INTENSITY  FUNCTION 


INTENSITY  FUNCTION 


« 


INTENSITY  FUNCTION 


INTENSITY  FUNCTION 


« 


4 


* 


214 


INTENSITY  FUNCTION 


L  INTRODUCTION 


The  scattering  of  light  by  aerosols  poses  interesting  theoretical  as  well  as 
experimental  problems.  The  atmosphere  contains  large  quantities  of  particular 
matter,  which  herein  are  called  "dust,”  that  affect  direct  (solar)  and  diffuse 
radiation.  Scattering  and  absorption  of  light  occur  not  only  in  the  ultraviolet 
(UV)  and  visible  spectral  range,  but  also  in  the  atmospheric  windows  of  the 
infrared  (IR).  Dust  consists  of  particles  which  differ  in  size  and  chemical 
composition.  The  natural  aerosol  population  in  the  troposphere  is  composed  of 
particles  originating  from  the  ground.  In  the  stratosphere,  Junge  (1961)  has 
suggested  that  the  natural  population  consists  of  three  families:  (a)  one  of 
tropospheric  origin  in  the  size  range  of  less  than  0.  1  micron,  (b)  one  of  strato¬ 
spheric  origin  in  the  size  range  between  0.  1  and  1.  0  micron,  and  (c)  one  of 
extraterrestrial  origin  in  the  size  range  greater  than  1.  0  micron.  As  to  the 
chemical  composition  of  stratospheric  aerosols,  particles  composed  of  ammonium 
persulfate  and  ammonium  sulfate  make  up  almost  all  of  the  aerosol  in  the  radius 
range  from  0.  1  to  1.  5  microns, which  particles,  in  turn,  make  up  more  than 
90  percent  of  the  total  mass  of  the  aerosol. 

Light  is  scattered  in  the  atmosphere  by  the  molecules  and  the  aerosol,  and 
attempts  have  been  made  to  infer  the  size  distribution  and  altitude  distribution 
of  aerosols  by  using  this  effect;  one  of  the  more  recent  is  that  of  Newkirk  and 
Eddy  (1963). 

To  interpret  experimental  results  based  on  scattering  of  light,  the  theoretical 
data  have  to  be  available  in  a  useful  form.  Nowadays,  scattering  functions  can 
be  computed  for  any  size  and  any  refractive  index  by  electronic  computers. 

Hence,  no  difficulty  arises  in  getting  the  basic  numerical  data;  but  this  is  only 
the  beginning  because,  in  designing  the  experiments,  it  is  helpful  to  have  a 
feeling  for  how  large  the  effect  of  aerosol  scattering  will  be.  In  some  cases, 
one  desires  to  maximize  the  scattering  effect;  in  other  cases,  one  wants  to 
minimize  it.  Furthermore,  it  is  helpful  to  know  the  angular  distribution  of  the 
scattering  intensity  and  polarization.  A  point  one  has  to  watch  out  for  is 
the  uniqueness  of  the  solution.  This  important  point  has  to  be  checked  very 
carefully.  Identical  scattering  coefficients  exist  for  different  size  parameters 
(within  particular  size  ranges)  because  of  the  oscillating  nature  of  the  functions;  ' 
or,to  express  it  in  physical  terms,  interference. 

The  purpose  of  this  contract  was  to  investigate  theoretical  scattering  data 
computed  by  electronic  computers  for  various  refractive  indices  and  size 
ranges  to  find  general  trends,  simple  but  for  experimental  work  accurate 
enough  interpolation  methods,  and  approximations  useful  for  rough  calculations. 
The  work  accomplished  has  been  reported  in  scientific  reports,  which  are 
summarized  in  this  final  report.  Revisions  and  additions  to  each  report  are 
given  if  new  information  was  obtained  after  the  report  was  completed,  or  if 
later  studies  necessitated  revisions  of  the  original  findings. 


TABU  IV  (cont. ) 


Definition  given  in  Table  II  or  III  not  used,  but  similar  parameter 
listed,  such  as  i^/of^,  K/of 

Huober  of  significant  digits  listed  (in  this  case  4,  e.g.  3.456  10* ^ ) 
Biufcer  of  places  after  decimal  point  (in  this  case  5,  e.g.  0.00023) 
Footnote,  appearing  at  end  of  table 


Refractive 

Ho. 

Index 

Range  of  06  Values 

Tabulated  Quantities 

\SS5S3S 

MEGGHjKHM 

36 

Boll 

o.a 

1-110  Ll8] 

K 

T,  G 

rel;  C, 

Guaprecht 

0.9 

1  -  .200  [20] 

K 

EHIAC;  4s 

Sllepcevlch 

0.93 

1  -  200  [25] 

K 

1954 

Oo 

1-90  [27] 

K 

37 

Gucker 

1954 

1.33 

3-3-18.5  [53] 

a,b,K 

•,  T 

rel. 

38 

Bailor 

1.05,1.1, 

0.2(0. 2 )7 

Tt§  (spec,  turbidity) 

T 

rel;  D,  5e 

Pangonls 

1.15,1.2, 

1954 

1.25,1-3 

39 

Johnson 

Terrell 

1.29 

,1  -  19.3  [22] 
0.5-25  W 

K 

T 

rel;  D,  3e 

1955 

K 

T 

40 

John* on 

1.33, 1.44 

0.5-6  [15]  # 

R  0-180 

T,  0 

rel;  3a, 

1955 

1.55 

based  on 

2.0 

R 

(23) 

4i 

Barker 

2.0 

3.3-12.5  til] 

1  0-40 

T 

rel;  D,  6s 

1955 

3. 3-4. 6  [5] 

a,b 

rel;  D,  6p 

1955U(L.C.) 

3.3-12.5 

a,b,i 

T 

*,  D 

42 

Penndorf 

1.0 

3^0.2(0.2)15, 

T 

rel;  D,  4p 

1956 

16-39-22  U3] 

K 

D,  4p 

2.0 

1.3-12.5  P3] 

K 

T 

— 

0.3-3. 2  W 

K 

T 

D,  4p 

43 

Penndorf 

1956 

1.33 

0.1(0.1)8(0.2)30 

(1)45 

K  (saootbed) 

T 

rel;  C  IBM 
701;  3s  based 
on  (44) 

44 

Penndorf 

1.33,1.4, 

0.1(0.1)30 

a,b,K 

€'  T  ,0 

rel;  C,  IBM 
701  10a  (a,b) 

5»(K) 

Goldberg 

1956 

1.44,1.486 

1.5 

45 

Bey,  Stewart, 

o«? 

0.0(0,01)10,# 

S.for  R  0-180 

•  T 

rel;  C 

Pinson,  Price 
1956 

*  # 

0 

5p 

46 

Lotnm 

1956,  U 

— 

1.1(0.6)9,5,10,20 

R  0-180 

■ 

47 

Siegel 

_ 

0.32,0.532,1.06, 

S,  R  0-0  -  180 

I 

1956,  U 

1.60,2.36,3.94, 

3.20,4.78 

i 

Author(() 
Taur  of  Wbl. 


KSSSSdm 

Index 


BSE» - 

Preaentetloo 


¥yp»  of  ~<ia*s>\ 
Accuracy 


of*  Value » 


tubulated  Qnentltlee 


77 


76 

79 


60 


81 


llengr-Bnttimn 

1962 


Rowell 

1962 

Delraendjian 

1963 


Puimdorf 

1962 

1963 


1.23 


0.1(0.1)10(1)20 

(3)50 


1.466 


1.29, 
1.315, 
1.525 

1  kV 

l.54,1.55| 
1.56 
2.2 


1-5 

1-33 


16  to  24 


0.5(0.5)15 

0.5(0.5)15 

0.5(0.5)15 

0.5(0.5)7 

0.5(0.5)10 

0.5(0.5)10 

0.5(0.5)10,12(4) 

40 

0-15 


0(0.1130 


1  0-0(2)10(10) 
170(2)180 


poaltlon  of  uuz  and 
■In  M>  to  160 


•  0-0,  160 


1,1,0-0-160 

1  0  -90 


C,  Bull 
dmmm  et 
ordlnutor 
4u,  5* 


ItMd  on 

(63) 

C,  1194  7090 
4a 


ST” 

C.Facom 
(126, 8» 


55 


3.2  List  of  Computations  arranged  according  to  Refractive  Index 


This  table  has  been  prepared  to  list  the  computations  for  increasing  re¬ 
fractive  indexes  only;  the  author  number  for  which  data  exist  has  been  in¬ 
cluded.  Those  being  considered  not  reliable  have  been  excluded. 

In  addition,  a  figure  shows  the  range  of  o  values  for  all  n  for  which 
data  exist.  Naturally  it  gives  only  a  general  picture,  but  even  that  might  be 
helpful. 
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with 


Zj  =  (n2  +  >^)  2  +  4  (n2  -X2)  +  4,  and 

Z2  =  4(n2  +  X2)2  +  12  (n2  -  X2)  +  9. 

The  total  scattering  coefficient  is  given  by 

=  - j[(n2  +  X2)2  +  n2-X2-2]  +  36  n"  *  , 

RR  3  z2 1  ’ 


2  X2  i  O4 


L,  6  f 

/n2  +  x2^2  -  4I  a2-  n  x  a3  1 

i  5  Zj 

in+xj  -e  a-  11  *  “  1 

V  J  3Zj  ' 

and  the  total  absorption  coefficient  as 
K  W  -  K_(e)  K>(8) 
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m.  CONTENT  OF  SCIENTIFIC  REPORTS 


Scientific  Report  1 ,  Results  of  an  Approximation  Method  to  the  Mie  Theory 
for  Colloidal  Spheres,  by  R.  Penndorf,  AFCRC-TN-59-608,  RAD-TR-59-38 
(September  1959),  43p, 

"An  approximation  method  is  outlined  for  computing  the  total  Mie  scatter¬ 
ing  coefficient  K  in  aerosols  using  infrared  light  sources.  The  method  is 
valid  for  any  size  parameter  o t,  from  small  to  very  large  spheres,  and 
any  arbitrary  real  refractive  index  n  <  2.  Analytical  expressions  are  de¬ 
veloped  for  the  phase  and  the  amplitude  of  K  at  the  maxima  and  minima. 

The  constants  in  these  analytical  expressions  are  determined  from  existing 
data  computed  by  the  exact  Mie  theory.  Finally,  a  graphical  interpolation 
process  allows  for  determination  of  K  for  any  arbitrary  size.  This  method 
leads  to  an  accuracy  of  ±2  or  3  percent  for  K  as  compared  with  results  ob¬ 
tained  by  the  Mie  theory. 

"This  approximation  method  is  than  applied  to  14  refractive  indices  n  between 
1.05  and  2.0.  The  results  for  smoothed  K  values  are  listed  for  the  range 
a  =  0.2  (0.2)  30  (0.5)  40,  which  is  considered  sufficient  for  most  applica¬ 
tions.  For  very  large  spheres,  the  laws  of  geometrical  optics  can  be  used 
successfully.  Graphical  representations  of  all  cases  are  also  included. 

From  these  figures,  the  appropriate  K  value  can  be  read  off  directly  for 
any  given  size  parameter,  and  also  for  any  given  radius  of  the  aerosol 
particle,  if  the  wavelength  of  the  incoming  infrared  radiation  is  selected." 

Scientific  Report  2,  Mie  Scattering  in  the  Forward  Area,  by  R.  Penndorf, 
AFCRC-TN-60-2^5,  RAD-TR-60-10  (February  I960),  87p. 

"The  scattering  of  light  is  investigated  from  a  theoretical  point  of  view. 

After  explaining  the  basic  Mie  formulas,  several  useful  scattering  functions 
are  defined,  such  as  the  angular  Mie  scattering  coefficient  ig  and  the  Mie 
phase  functions  p^j.  Next,  the  scattering  pattern  for  the  forward^  sideward, 
and  backward  area  is  briefly  described,  followed  by  a  very  detailed  investi¬ 
gation  for  the  forward  direction,  which  is  based  on  all  available  data  com¬ 
puted  according  to  the  Mie  formulas.  A  wide  range  of  refractive  indices  n  . 
from  n  =  1.05  to  n  =  2.0  has  been  selected,  and  all  values  of  ig  and  p^j  are 
extensively  listed  in  the  appendix.  The  behavior  of  these  scattering  coeffi¬ 
cients  is  investigated  and  described  in  detail.  It  is  found  that,  for  example, 
while  ig  shows  major  oscillations,  only  plateaus  are  found  for  pM,  and  their 
width  is  related  to  the  refractive  index.  The  amplitude  of  the  major  oscilla¬ 
tion  decreases  with  increasing  size  parameter.  Superimposed  on  these  os¬ 
cillations  are  ripples,  the  phase  of  which  agrees  exactly  with  the  phase  of 
the  ripples  for  the  total  scattering  coefficient  K.  Several  ripple  systems 
are  noted,  and  some  regularity  in  their  appearance  is  clearly  indicated. 

The  multi  value  ness  of  the  scattering  function  and  its  consequences  are  also 
pointed  out. 
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Corrections 


page  15,  paragraph  3,  the  statement  made  is  wrong,  van  de  Hulst's 
formula  for  is  correct.  It  leads  to  the  correct  formula  in  agree¬ 
ment  with  equation  (17a),  provided  that  n  =  n  -  ix  is  inserted  for  m  in 
van  de  Hulst's  formula. 

page  44,  paragraph  4,  last  sentence,  change  signs  to  read:  "For  the 
chosen  values,  the  error  ranges  for  a=  1  from  -11  to  -30  percent 
for  n  =  1.  25,  from  3  to  -6  percent  for  n  =  1.  5,  and  from  10  to  16  per¬ 
cent  for  n  =  1.  75.  " 

page  46,  Figure  11,  the  three  refractive  indices  are  n  =  1.  25,  1.  50, 
and  1.  75. 

Scientific  Report  4,  Bright  and  Dark  Rings,  by  R.  Penndorf,  AFCRL-425,  RAD- 
TR-61-16  (March  1961),  24p. 

"The  scattering  diagrams  of  spherical  aerosols  possess  maxima  and  minima, 
which  can  be  interpreted  as  sequences  of  bright  and  dark  rings.  These 
rings  are  investigated  for  n  =  1.10,  1.20,  1.30,  1.33,  1.40,  1.44,  and  1.50. 
The  results  show  a  systematic  behavior  exists  for  the  appearance  of  these 
rings,  such  as  the  number  of  rings  and  their  regular  shift  in  angular  posi¬ 
tion  as  function  of  the  size  parameter  a.  For  the  intensity  function  ij,  the 
rings  due  to  diffraction  move  with  increasing  size  parameter  into  the  for¬ 
ward  scattering  area,  whereas  those  caused  by  refraction  and  reflection 
within  the  sphere  stay  in  the  backward  scattering  area.  The  rings  appear 
first  in  the  backward  area  as  long  as  the  size  parameter  is  small,  but  from 
p  >  4  (p  =  normalized  size  parameter)  on  they  begin  to  appear  at  80  to  90 
degrees;  i.e.,  in  the  right-angle  scattering  area.  For  the  intensity  func¬ 
tion  i2  however,  the  systematic  behavior  is  poor  and  ill  defined.  Never¬ 
theless,  some  of  the  basic  features  are  similar  to  those  for  ij .  Finally, 
a  system  similar  to  ij  exists  for  the  angular  Mie  scattering  coefficient  ig. 

In  the  last  case,  data  are  given  for  n  =  1.33  only." 

Corrections 

page  5,  Figure  5,  the  minima  possess  no  cusp  as  indicated.  The  mini¬ 
ma  should  show  a  smooth  transition  occurring  somewhat  higher  than 
indicated.  The  same  is  true  for  figures  3  and  4. 

page  10,  Figure  5-7  should  be  revised,  see  the  discussion  in  Scientific 
Report  10,  page  7.  We  believe  now  that  our  first  interpretation,  as 
given  in  this  report,  is  not  always  correct.  There  are  minima  which 
"vanish"  as  a  increases.  The  better  solution  is  indicated  by  figures  3-6 
and  3-7  in  Scientific  Report  10,  and  in  section  II  of  this  final  report. 


-23- 


ing  angle  in  the  forward  direction  ie  defined  ae  0  -  0°,  in  agreement  with 
v.d.  Bulat  and  the  present-day  use.  In  the  older  literature  it  will  be 
found  that  O'  “  l80°  is  defined  as  the  forward  direction,  which  goes  back 
to  Mie's  definition.  However,  the  use  of  0  =  0°  for  forward  scattering  is 
advisable,  and  our  tables  are  listing  the  scattering  angle  in  our  defini¬ 
tion  irrespective  of  what  definition  the  author  used. 

Since  this  mix-up  introduces  difficulties  in  using  the  literature 
properly,  a  special  symbol  (0)  is  added  in  our  lists  to  indicate  to  the 
reader  that  he  has  to  look  out  for  the  definitions  used  by  the  author. 

Some  authors  do  not  tabulate  i, ,  ig  directly,  but  related  quanti¬ 
ties  such  as  i^o^,  ie  ■  (l^+ig )/27f'c^,  phase  functions,  or  other  quantities. 
In  those  cases  we  have  used  the  notation  i^,  ig  and  qualified  it  by  the 
syrtol  (#)  to  indicate  to  the  reader  that  a  related  quantity  is  tabulated. 

4.  The  me  functions  are  a^,  b^.  They  sure  useful  if  additional 
computations  are  required,  problems  arise  for  tvo  reasons,  first  the 
definition  of  forward  scattering  o'  «  180°  instead  of  0  ■*  0°  introduces  a 

constant  factor.  If  a  is  the  me  function  using  8*0°  for  forward  scat- 

n 

taring  a^  that  for  d'  -  180°  for  forward  scattering,  both  are  related 


by 


(-1)“  1 


(-D-1  i  b* 

m 


<*0 


2m  +  1 


This  change  ie  Important.  Secondly,  the  term  y  bas  been  included 
In  the  definition  of  a  ,  for  example  by  Lowan  and  by  Guaprecht  and  SUepcevlch. 

Benoe,  their  a  and  b  must  deviate  from  those  given  in  Table  II  (Bq.  2.5). 

b  n 

5.  The  radar  cross-section  cr .  The  tables  indicate  by  the  letter 


B  that  this  quantity  has  been  computed.  It  is  an  unfortunate  definition  but 
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IV.  ABBORBHIO  BTHBRE8 


4.1  List  of  existing  computations  for  »m,  b^,  8.^  S2,  i1#  i2,  K*,  K*,  K* 

and  related  functions  or  coefficients  for  nonodliperse  Isotropic 
spherical  particles  (dielectric  parti clea). 

The  following  table  requires  again  soma  explanation.  The  first  column 
lists  the  running  number.  This  number  agrees  with  the  nuober  in  the  refer- 
ence,  section  4.3*  The  second  coluian  lists  the  "name  of  the  author (■ )  and 
the  year  of  publication".  The  symbol  U  after  the  year  means  that  the  data 
are  unpublished  or  deposited  as  document  in  the  Library  of  Congress.  The 
third  coluian  contains  the  refractive  index *nj  according  to  Bq.  (1).  Values 
froa  authors  who  used  Bq.  (2)  and  (3)  have  been  transformed  to  our  system. 

Is  some  cases,  if  it  has  not  been  stated  accurately,  no  data  are  given  in 
this  column. 

The  fourth  coluan  lists  the  type  of  compound  to  which  it  applies,  such 
as  Au  la  HgO,  or  carbon,  or  water  droplets  in  air  (HgO). 

The  fifth  columa  shows  the  appropriate  spectral  range  for  which  this 
value  of T and  the  compound  apply.  V  (visible),  Ut,  (infrared),  R  (radar) 
are  the  sy*ols  used.  They  seen  to  be  clear.  Sometimes  the  wavelength  is 
also  listed  for  V  and  IR  the  numbers  are  given  la^,  for  R  in  cm. 

The  sixth  column  shows  the  "range  of  (X  values".  This  system  is  the 
same  as  in  section  3.1.  The  seventh  colum  shows  the  "tabulated  quantities". 
The  eighth  and  ninth  columi  is  again  similar  to  those  in  section  3.1. 

Data  for  some  water  and  ice  mixtures  and  pure  water  and  ice  are 
given  in  section  VI. 
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6.2  ATLAS  OF  SCATTERING  DIAGRAMS  FOR  n  =  1.1,  1.2,  1.33,  1.4,  1.44,  and  1.5 


Bie  diagrams  for  n  -  1.1  and  n  =  1.2  are  based  on  numerical  data  published 
by  Pangonis  and  Heller,  the  rest  is  taken  from  our  unpublished  listings.  The 
diagrams  for  n  =  1.33  and  n  =  1.5  are  repeated  from  earlier  reports. 
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